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ABSTRACT 

We present spatially- resolved K- and L-band spectra (at spectral resolution 
R = 230 and R = 60, respectively) of MWC 419, a Herbig Ae/Be star. The 
data were obtained simultaneously with a new configuration of the 85-m baseline 
Keck Interferometer. Our observations are sensitive to the radial distribution of 
temperature in the inner region of the disk of MWC 419. We fit the visibility 
data with both simple geometric and more physical disk models. The geometric 
models (uniform disk and Gaussian) show that the apparent size increases linearly 
with wavelength in the 2-4 /im wavelength region, suggesting that the disk is 
extended with a temperature gradient. A model having a power- law temperature 
gradient with radius simultaneously fits our interferometric measurements and 
the spectral energy distribution data from the literature. The slope of the power- 
law is close to that expected from an optically thick disk. Our spectrally dispersed 
interferometric measurements include the Br 7 emission line. The measured disk 
size at and around Br 7 suggests that emitting hydrogen gas is located inside (or 
within the inner regions) of the dust disk. 

Keywords: stars: individual (MWC 419); stars: pre-main sequence; (stars:) 
circumstellar matter; stars: emission-line, Be; techniques: interferometric; instrumentation: 
interferometers 

1. Introduction 

Pre-main sequence (PMS) stars fall to the upper right of the main sequence on the 
Hertzsprung-Russell (HR) diagram, journeying towards the main sequence via radial 
contraction. For the first several million years, a PMS star is surrounded by a disk of gas 
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and dust, left over from the early stage of star formation. 

The evolution of a disk is not well understood. Circumstellar disks provide the 
raw material for planet formation. Thus, understanding the evolution of a disk helps us 
understand planet formation. Clues to the physical conditions of planet formation and 
hence future suitability for life on planets other than Earth potentially can be inferred from 
a detailed characterization of inner young stellar object (YSO) disks. 

Herbig Ae/Be (HAeBe) stars are intermediate- mass PMS stars with spectral types 
earlier than GO, located close to the zero-age main-sequence . Early spectral energy 



distribution (SED) modeling efforts ([Hillenbrand et al. 1992T ) to explain the infrared 



excess of HAeBe stars assumed (a) optically thick but geometrically thin circumstellar 
disks with inner optically thin holes, up to several stellar radii, in order to account for 
the observed inflections in the 1-5/zm region, and (b) relatively large accretion rates, on 
the order of 1(T 6 M yr 1 , in order to account for the str ength of the 3.8/im (L-band) 



excess. It was quickly pointed out (IHartmann et al. 19931 ) that at such a large mass 



accretion rate, the gas in the inner dust-free hole is expected to be optically-thick and 
to produce excess near-infrared emission, which is inconsistent with the dip seen in 
the near- infrared SED. For the relatively low accretion rates later inferred for HAeBe 
objects, disk models could not explain the observe d SEDs in the nea r- infrared region. 



Puffed-up inner dust-rim models were intro duced (INatta et al. 20011 ) and later refined 



flDullemond et al. 2001 



Isella &: Natta 20051 ) to ameliorate the shortcomings of the classical 



disk models. These later models attribute the near-infrared excess to stellar radiation 
shining directly on the inner dust edge of the dust disk. Tuthill et al. (2001) directly 
imaged the high luminosity YSO LkHa 101 and independently proposed that the bulk of 
the near-infrared emission arises in a hot ring located at the dust sublimation radius. 

Long-baseline interferometric observations provide the milliarcsecond angular 
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resolution required to resolve the planet-forming structures immediately surrounding 
PMS stars. Such high-angular-resolution infrared observations are well suited for probing 
near-circumstellar environments since the inner zones of the circumstellar disks and 
envelopes emit primarily at near-infrared wavelengths. A large number of YSOs across 
the luminosity range have now been spatially r esolved at near-in f rared wavelengths 



(IMillan-Gabet et al. 2001 



Colavita et al. 20031: Eisner et al 20031: Eisner et al. 2004 : 



Eisner et al. 2007 



Monnier et al. 2005: Akeson et al. 2005: Millan-Gabet et al. 2007 



Isella et al. 2008 



Eisner et al. 2009 



Tannirkulam et al. 20081 ). showing that indeed the 



characteristic sizes correlate strongly with central luminosity, lending support to the 



puffed-up" inner dust rim paradigm, especially 



also pointed out (IMonnier fc Millan-Gabet 2002 



or Herbig Ae objects. However, it was 



Eisner et al. 2004 



Monnier et al. 20051 ) 



that most higher luminosity Herbig Be objects are considerably undersized compared to 
the predictions of this model, and in better agreement with the "classical" models featuring 
optically thick emission that extends inward very close to the central star. 

MWC 419 is a photometrically variable emission line star th at is in many ways typical 



has spectral typ e B8 ( 



nebula ( 


Herbig 1960 




V1WC 419 
), and 


(Hillenbrand et al. 1992 



distance 650 pc (jHillenbrand et al. 19921 ) ; it is also known as V594 Cas or BD+61 154. 
Narrow band Ha imaging st udies show unipolar large scale structures around MWC 419 
(IMarston fc McCollum 20081 ). Interestingly, these images show two lobes on the south-east 
side of MWC 419 - the outer lobe extending to ~ 4 arcminutes - suggesting a history of 
episodic mass-loss in this PMS star's evolution. These authors attribute the non-detection 
of north-west lobes to obscuration. These large scale structures have no influence on 
our observations of the c entral region. The mea sured intrinsic polarization of MWC 419 



in the V-band is 0.53% (jHillenbrand et al. 1992l ) suggesting that the contributions from 



reflected optical light to the total emission is insignificant. No close companions to 



- 5 - 



MWC 419 are known flPirzkal et al. 19971 : Bernacca et al. 19931 ). The star has P-Cygni 
characteristics in the lower Balmer lines and the implied wind has been modeled in detail 
by Bouret & Catala (1998). However, it has also been considered as a B[e] star and 
interpreted as undergoing post-main sequence evolution rather than pre-main sequence 
evolution. T 



r e star is projected on the outer regions of the young association NGC 225 



(IHerbig 1960 



CO emi ssion ( 



Subra maniam et al. 20061). within a small area of enhanced extinction and 



Hillenbrand 1995 



Fuente et al. 20021 ) that is associated with the dark cloud 



L 1302 (Ivan den Ancker et al. 19981 ). However, the membership of the star in the cluster 
i s debated (Lattanzi et al. 1991). The age of the cluster is often quoted as 120 Myr 
( jLattanzi et al. 199ll ) but Subramaniam et al. (2006) have argued based on the presence of 
a number of young stars that it is really a <10 Myr pre-main sequence population. 

As an IRAS and ISO source, the infrared dust spectrum an d spectral energy 



distribution of MWC 419 have been discussed by several au t 



Chen et al. 2000: 


Berrilli et al. 1992: 


Pezzuto et al. 19971). No Spitzer dat 



Hillenbrand et al. 1992 



rors ( 



Brooke et al. 1993 



Miroshnichenko et al. 1997 : 



silicate emission at 10 /zm in ground-based data (IChen et al. 20001 ) . 



MWC 419 was previously observed interferometrically using the Palomar Testbed 
Interferometer (PTI), and was spatially resolved in the K-band with a reported uniform-disk 
angular diamete r of 3.34 ±0.16 milliarc sec (mas), or a Gaussian distribution FWHM of 
2.07±0.11 mas flWilkin & Akeson 20031 ). MWC 419 was not resolved from the Infrared 



Optical Telescope Array (IOTA) interferometer observations (IMillan-Gabet et al. 20011 ) . 
where the angular resolution was ~ 6 mas at 2.2/im. 

The Keck Interferometer (KI), with ~3 mas angular resolution at 2.2/im, can resolve 
the inner disks (~ 1 AU) of nearby PMS stars. In this paper we report observations of 
MWC 419 using multi-color interferometry at well-separated wavelengths, here K-band 



- 6- 



(2.0-2.4 //m) and L-band (3.5-4.1 /im). The simultaneous K- and L-band interferometric 
measurements enable us to probe different regions (~ 1300K and 800K respectively) of the 
inner circumstellar disks of PMS stars. Discrete spatial distributions, such as dust-rims, 
and relatively-smooth spatial distributions, such as classical accretion disks, are expected 
to have different size vs. wavelength behaviors and can be distinguished in such multi-color 
observations at well separated wavelengths. In addition, interferometric measurements in 
the relatively unexplored L-band provide further constraints to the disk/envelope geometry 
via temperature- and density-sensitivity of different models. Multicolor interferometry at 
well separated wavelength s (H, K and N-ba nds) was recentl y reported for two other YSO 



disks, namely, MWC 297 flAcke et al. 20081 ) and MWC 147 flKraus et al. 20081 ) from VLTI 
observations. 

In Section 2, below, we present our observations and data reduction. In Section 3 
we describe the data analysis where we fit the visibility-squared (V 2 ) data with various 
YSO disk models. In Section 4 we discuss our results and in Section 5 we provide a brief 
summary. 



2. Observations and Data Reduction 

KI is a near- and mid-infrared long-baseline interferometer consisting of two 10-m 
diameter apertures separated by a B = 85-m baseline at a position angle of ~ 38° east of 
north. Both Keck telescopes are equipped with adaptive optics. The maximum resolution of 
our KI observations is X/2B ~2.7 mas and ~4.5 mas in the K-band and L-band respectively. 

The L-band science instrument, consisting of an infrared camera and beam-combiner 
optics, is similar to the existi ng H- and K-band in strument (FATCAT) routinely used for 



V 2 measurements since 2002 fjColavita et al. 20031 ). We obtained first light with this new 
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instrument in April 2008. The L-band limiting magnitude of 6 is set by the requirement 
for a broadband phase signal-to-noise ratio (SNR) of 10 in the 10 ms fringe-tracker frame. 
The total integration time on source is typically 200 sec, which provides good SNR for V 2 
measurements in the 10 spectral channels. 

The L-band instrument uses a PICNIC focal plane array detector with a 5 /zm cutoff 
wavelength, while FATCAT uses a HAWAII array. The detector electronics are similar for 
both instruments. Details of this L-band instrument are given in Ragland et al. (2008). 

For the current work a prototype multi-wavelength observing capability was used, 
enabling simultaneous K- and L-band observations of our science target. In this 
configuration, the telescope pupil is split into left and right-halves at the dual-star modules 
of both telescopes and routed through separate coude paths utilizing the existing beam-train 
infrastructure for the nulling mode. The left-half pupils of both telescopes are combined by 
the K-band instrument and the right-half pupils are combined by L-band instrument. 

The K- and L-band science instruments each have 2 complementary interferometric 
outputs with different spectrometers on each. For L-band, these are a (pseudo) broad-band, 
and a low dispersed mode (10 channels across the L-band; R = 60). For K band, these are 
a broad-band, and a medium dispersion mode (42 channels across the K-band; R = 230). 

The field-of-view of the instrument, defined by the single mode fibers that couple that 
light to the detector array, is ~ 55 mas for the K-band and ~ 93 mas for the L-band 
observations. These field restrictions were enforced in our modeling work. 



The ZABCD algorithm (jColavita et al. 20031 ) is used for fringe tracking and science 



measurements. In this procedure five reads are made while the fast delay line scans over one 
wavelength (2.2 /zm for the K-band and 3.7 /xm for the L-band); explicitly, the reads are the 
reset pedestal (z), followed by four non-destructive reads (a, b, c, d) spaced at A/4 intervals 
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as the detector integrates up. For each quarter-wave bin we then calculate A = a — z, 
B = b — a, etc. These values are used to estimate the square of the fringe visibility (V 2 ) 
for science, and the fringe phase and SNR for fringe tracking. We apply corrections to the 
bin data to account for detector bias and for differences between the length of the delay 
scan and the actual wavelength. Residual instrument bias and atmospheric seeing effects 
are corrected by observing a calibrator of known V 2 under similar observing conditions. 
Essentially, we estimate instrument transfer function using the bracketing observations of 
calibrator stars with known angular diameters using a weighted averaging scheme that 
considers the time and sky proximity of calibrators relative to the target 1 . 

The observations reported here were taken on the night of UT 19 August 2008. We 
observed two calibrators - HD 1843 and HD 6210 - under similar observing conditions as 
the science target to calibrate the science dat a. The adopted angular sizes of the calibrators 



are 1.0±0.02 and 0.5±0.01 mas respectively (Ivan Belle 19991 ). We performed bracketed 
calibration meaning that our observing sequence consisted of calibrator-target-calibrator 
measurements. 

The measurements presented here were taken over a narrow range of position angles 
(25° to 33°) and projected baselines (76.1m to 77.8m), and hence no attempt was made to 
derive an inclination angle and position angle of the disk. For the purpose of this article, 
we assume that the disk is face-on. Thus, the sizes reported here are the values along the 
position angle of about 29°, and the actual size of the disk could be larger depending on 
the inclination angle and position angle. Our K-band measurements, taken at the position 
angle of about 29°, gives similar disk size as that of previous PTI observations, taken at the 
position angle of about 83°. Hence, a face-on disk for MWC 419 is a resonable assumption 
in the absence of necessary high spatial resolution observations. 
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3. Analysis 

In this section, we fit our measurements with different models of increasing complexity 
and increasing physical realism. 

The measured V 2 includes contributions from the central star (V*) and the bright 
circumstellar disk. The visibility-square (V^ sfc ) of the circumstellar disk of MWC 419 is 
obtained from the measur ed data by removing the contributions from the central star 



flMillan-Gabet et al. 20011 ) using 



/ FjfVtf -\- FdiskVdisk \ \ 

' measured I Z7 i E" 1 / \ / 

\ r ^ + fdisk J 

The adopted disk-to-star flux ratio (F^sfc/F*) is 12 and 40 in K and L based on our 
SED analysis discussed in Section 3.3. The central star is assumed to be unresolved for our 
observations (i.e., V* = 1.0), which is a reasonable assumption for a B8 star at a distance 
of 650pc. These corrections yield values of V^ isk that are smaller than the total V 2 by 
0.04 (~ 9%) and 0.01 (~ 2%) in K and L respectively. The K-band V 2 measurements 
are consistent with earlier K broadband measurements and provide additional spectrally 
resolved information within the K-band spectral region. 

The visibilities of the disk models (Section 3.2 and 3.3) are computed by numerically 
summing the contributions from annular rings of infinitesimally small widths and weighting 
them by their respective flux contributions; for this purpose, we divided the disk radially 
into 5000 annular rings, in logarithmic scale. The normalized visibility for a uniformly bright 
annular ring with an inner diameter of 9 in , an outer diameter of 6 out and in < 9 < 9 out , can 
be written as 



V ann {B, A, 6) = I „ //l2 2A „ 2 . ) (OavtJiinBOart/X) - ^(ttS^/A)) (2) 



nB(d 2 out -e 2 n ) 
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where J\ is the first-order Bessel function, B is the projected baseline and A is the 
wavelength of observation. 

The flux from such an annular uniform ring can be written as 



F ann (\,0) = ^P(T,\) (C-O (3) 

where P(T, A) is the Planck blackbody function and T is the mean temperature of the 
annular ring. 



Now, the visibility of the total disk can be written as 



Vdisk(B, A) = ^2 f Krni(-B) \ 



6) x _f°™( A ' g ) ) ( ii 



Fann ( A , 

where F ann (X, 9) is the total disk flux derived by summing all the annular rings of the 
disk in order to normalize the flux. We adopted a nonlinear least-squares fitting method to 
minimize the chi-square for all our modeling work. 

For our SED analysis, photometric measurements from the literature are used, after 
interstellar extinction corrections using the extinction law of Car delli et al. (1989). The 



extinction for MWC 419 in the V-band is assumed to be A v = 2.1 ([Hillenbrand et al. 19921 ). 



The models presented in this article assume face-on geometry for the disk. A disk with 
an inclination angle <fi could predict a larger inner disk radius depending on the inclination 
angle and the position angle of the projected interferometer baseline on the sky. In this 
case, SED models would underestimate the flux at all wavelengths by a factor of cos(<f>). In 
other words, classical accretion disk models would have a factor of cos(0) larger accretion 
rate to explain the SED data and power-law models would have a relatively larger disk 
temperature at all radial distances. 
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3.1. Simple wavelength-dependent geometrical models 



The measurements presented in this article consist of spectrally dispersed data within 
the K- and L-bands. As stated above, the K and L visibilities corrected for the stellar 
contribution are different, suggesting a wavelength dependence to the size of the spatially 
resolved emission. We can use the spectrally dispersed data to investigate in more detail 
this wavelength dependency. We choose three geometrical models, namely, uniform-disk, 
Gaussian distribution and ring models, to fit to our measurements. 

The complex visibilities for a normalized pole-on uniform-disk and a Gaussian 
distribution can be written as follows: 



Here Our, is the uniform-disk angular diameter of the circumstellar disk, 9 FWHM is the 
full-width at half maximum (FWHM) of the Gaussian distribution, and other terms as 
above. See Eqn. [2] for the ring model. 

To start, we fit the data with a wavelength-independent geometrical model (Figured]). 
The derived uniform-disk, Gaussian and ring angular sizes, from simultaneous fits to K 
and L band data, are 3.59 ± 0.01 mas, 2.21 ± 0.01 mas (FWHM) and 3.91 ± 0.01 mas 
respectively. The corresponding reduced-chi-square (xr) values are 36.1, 35.0 and 38.3 
respectively. As indicated by the very poor \r values, the uniform disk, Gaussian and ring 
models fit neither the K and L data in the mean, nor the slope of the visibilities throughout 
either the K- or L-bands. 

As the wavelength-independent size models fail to fit the measurements, we fit the 



V UD (B,X) 



2J 1 [nB9 UD /\] 
ttB6 ud /\ 



(5) 
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Wavelength Independent Models 

qqq r i 




20 25 30 35 40 
Spatial Frequency (10~ 6 B/X) 



Fig. 1. — Wavelength- independent model fits. The mult i- wavelength measurements (dia- 
mond symbols) are plotted against spatial frequency. Also shown are geometrical model fits 
to the data; the dotted-line, dashed-line and continuous line refer to uniform-disk, Gaussian 
and ring models respectively. The poor model fits indicate that these simple models are an 
inadequate description of the data. 
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data with wavelength-dependent sizes. The results of these geometrical model fits to the 
measured data points are shown in Figure [2J 



= 1.?07 + (firr 



0(rras) = 0.790 + 0.612A (firr. 



Fig. 2. — Wavelength-dependent linear fits. Left: The derived uniform-disk angular sizes as 
a function of wavelength are shown in diamond symbol along with error bars. A linear fit 
to these sizes (dotted line) is also shown here. The dip seen at around 2.17/im is due to the 
presence of compact Br 7 emission line (see Section 4 for details). Middle: Same as the left 
figure, but a zoomed view of the Br 7 emission line region. Right: Same as the left figure 
except for a Gaussian distribution. 

The observed wavelength dependence in the 2-4 /im region has a simple linear 
relationship to first order. The parameters of a linear fit to the derived apparent wavelength 
dependent diameters are given in Figure [2J The Xr value for this model fit is 2.1. The 
measured uniform-disk diameter in the center of the L-band (8l = 5.04 mas) is ~ 44% 
larger than in the center of the K-band {6k — 3.49 mas). The steep slopes of these linear 
relations suggest that the 2-4/im emission source must be extended with strong radial 
temperature dependence. Earlier multi-color (near- an d mid-infrared) interferometry of tw o 



Kraus et al. 2008^ . 



YSO disks have shown similar wavelength dependence (lAcke et al. 2 008; 
In the following section we explore more complex extended disk models to explain our 
observations. 
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3.2. Geometrically thin, optically thick accretion disk models 

We compared our interferometric data with an accretion disk model by Hillenbrand et 
al. (1992) based on an SED analysis. In this classical accretion disk model, the temperature 
distribution of dust is derived by combining the contributions from a radiatively heated 
reprocessing component T rep and a viscously heated accretion component T acc . The 
parameters of the disk are the following: The inner disk (hole) radius is Rj n = 0.22 AU (10 
R*), the outer disk radius is 62 AU and the mass accretion rate (M) is 1.98 x 10~ 5 M 
yr _1 . The stellar radius is 4.8R , the stellar effective temperature is 11,220K, the distance 
to the star is 650pc and the stellar mass is 5.3M . The inner disk temperature is 2480K 
and the outer disk temperature is 40K. The resulting flat accretion disk model highly 
overestimates V 2 in the K- and L-bands. We get a reduced-chi-square (xr) value of 664 for 
the interferometric data. The model results are given in Table 1. The physical reason that 
the classical accretion disk model fails is that it creates too much flux in the inner region of 
the disk because of the added accretion luminosity, making the angular size much smaller 
than observed. 

We have also fit visibility data with a similar accretion disk model (Figure [3D by 
treating the inner disk (hole) radius as a free model parameter. In addition, the mass 
accretion rate was increased by a factor of 1.5 (i.e. 2.97 x 10~ 5 M Q yr _1 ) in order 
to fit SED data. The \\ value for this model fit to interferometric measurements is 
2.59. The derived inner disk (hole) radius is 1.66 ± 0.01 mas (0.54AU; ~24R iiI ) and 
the inner disk temperature is 1457K. The outer disk radius is fixed at 50 AU where 
the disk temperature is 51K. This model fits the K band data well, but is about 
3 a too high in V 2 compared to the L band data. The extremely large accretion 
luminosity required by the SED data poses challenges to these two accretion disk 
models. Thus the classical model that could fit the SED fails here, just like in the previous 
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interferometric studies of Herbig Ae/Be disks (IMillan-Gabet et al. 200 It Eisner et al. 2004 ) 
Moreover, it is not sufficient to adjust the inner dust radius and temperature alone. A 
modification of the disk temperature profile is needed a s in the case of more recent studies 



(Eisner et al. 2007 



Acke et al. 2008 



Kraus et al. 20081 ); this is discussed in the following 



section. 
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Fig. 3. — Classical accretion disk model fits. Interferometric data points are shown with 
error bars. Shown in continuous line is an accretion disk model with the inner radius of the 
disk (hole size) treated as a free model parameter. The outer radius of the disk is fixed at 
50 AU. Right: Photometric data taken from the literature are shown along with the SED 
model for the same accretion disk model. The dashed line shows the blackbody spectral 
distribution of the central star and the dashed-dotted line shows the blackbody emission 
from the disk. 
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3.3. Geometrically thin, optically thick disk with power-law temperature 

distribution 

The wavelength dependence of the measured visibilities and SED data shown above 
implies an inverse T(r) relationship. We fit our data with a simple model with a power-law 
temperature gradient of functional form T(r) oc r~ a , where r is the radial distance from 
the central star and a is the power-law parameter. The model V 2 is derived using Eqn. [2] 
and Eqn. [3] with a power-law temperature distribution. The radius of the inner disk and a 
are treated as free parameters. The temperature of the inner disk is fixed at 1800K and 
the optical depth at 1 for the entire disk, in order to satisfy the SED data; a larger value 
for the optical depth overestimates the disk flux. The resultant model fits are shown in 
Figure HI The derived inner disk radius is 1.47 ± 0.02 mas (0.477AU; ~ 21R*) and the 
derived value for a is 0.71 ± 0.01. The outer radius is fixed at 20AU (in order for the SED 
to be consistent with the far-infrared photometric data) and the corresponding temperature 
is 125K. The \r f° r this model fit to our interferometric data is 0.73, and the sensitivity 
of Xr to the model parameters is shown in Figure [5j As the power-law disk model is a 
satisfactory fit, we consider it the simplest model that explains the data. However, one can 
also conceive of more complex and physically plausible models such as optically thin inner 
disk regions (holes) surrounded by optically thick outer disks, perhaps including puffed up 
inner rims. Our initial exploration of such models suggests that they can also fit our data, 
but the interferometric data themselves do not drive us to such complex models. 

3.4. Binary star model 

In this section we investigate whether the data can be explained, without invoking 
hot/warm circumstellar dust, by the effect of a companion star. The visibility model for a 
binary system can be written as 
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Fig. 4. — Power-law disk model fits. Left: Interferometric data points are shown with 
error bars. The solid line is a power-law temperature gradient disk model. The power-law 
parameter, the inner disk radius and the inner disk temperature are free parameters. The 
outer radius is fixed at 20AU. Right: Photometric data taken from the literature are shown 
with error bars. The solid line is the SED model for the same accretion disk model. The 
dashed and dashed-dotted lines are the SEDs of the star and disk. 
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Model 


Rjn 


T 

-L in 


a 


M 


x\ 




(AU) 


(K) 




(x 10~ 5 M yr" 1 ) 




Classical accretion disk 


0.22 


2480 


0.75 


1.98 


664 


(parameters fixed at published values) 












Classical accretion disk 


0.538±0.001 


1457 


0.75 


2.97 


2.59 


(Rj n and M varied) 












Modified power-law disk 


0.477±0.005 


1800 


0.71±0.01 




0.73 



Table 1: Derived parameters for the disks models presented in Sections 3.2 and 3.3. The \\ 
values are for the interferometer data (only). 



1 http : / /nexsci . caltech . edu / software/ V2calib / 

yjvf + R 2 V s 2 + 2 jig |K| R cos((2n/\) B ■ s) 

VbinaryW = 1+ R 

where R is secondary-to-primary flux ratio, V p and V s are visibilities of primary and 
secondary components, and s is the binary separation. 

The derived parameters from a binary model fit to our measurements are the following: 
projected binary separation along the position angle of 29° (east of north) = 5.12±0.01mas, 
primary-to-secondary flux ratio = 11.79±0.30, and uniform-disk angular diameter of the 
primary component (YSO disk size) = 3.53±0.01 mas (a binary model with two point 
sources fails to fit SED data). The diameter of the secondary component is fixed at 0.001 
mas (effectively unresolved). The model fit to the data is shown in Figure [6J The x\ value 
is 0.94. In this scenario, the secondary component would have the same K-L color of 1.34 
as the primary YSO disk suggesting that the potential companion were a YSO disk with 
significant IR excess. In this case, our observations would have resolved the disk around 
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\2 for Power— law models 

q 54 r 

0.52 - — - 




o.42 r 1 

0.65 0.70 0.75 0.80 

a parameter of the Power— law 

Fig. 5. — Contour map of Xr as a function of the power- law parameter a and the inner disk 
radius. The contour lines - from inner to the outer - refer to \\ values of 1.0, 1.2, 1.4, 1.6, 
1.8 and 2.0. 

the companion. Hence, we also fit our measurements by fixing the companion size to be 2 
mas in order to explore this possibility. The derived parameters for this binary model are 
the following: projected binary separation along the position angle of 29° (east of north) 
= 5.23±0.01mas, primary-to-secondary flux ratio = 11.6±0.3, and uniform-disk angular 
diameter of the primary component (YSO disk size) = 3.51±0.01 mas. The model fit to the 
data is also shown in Figure El The Xr value is 1.57, but additional observations at other 
baseline orientations would provide significant additional constraints on the binary model. 
The derived flux ratio of ~ 12 suggests a cooler stellar companion with a smaller inner disk 
radius compared to the primary star. However, such a companion model would overestimate 
SED at the shorter wavelengths because of the flux from the companion photosphere. 
Hence, a binary scenario is not favored as an explanation for our measurements. 
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Fig. 6. — Binary model fits. Left: The multi- wavelength measurements are plotted against 
spatial frequency. Also shown is the best fit binary model (continuous line) assuming an 
unresolved companion. Right: Same as the previous caption but for a companion of 2 mas 
diameter. 
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4. Discussion 



The KI measurements pre sented here agree with p revious observations from PTI of 



the K-band size of MWC 419 flWilkin fc Akeson 20031 ). The addition of the spectrally 



dispersed K-band information and new L-band observations, also spectrally dispersed, 
provide powerful new constraints on the physical structure of the material surrounding the 
central star. As shown the left-most panel of Figure 2, the wavelength dependence of the 
size is apparent within both the K and L bands; the combination of the two bands allows 
an even tighter constraint on the temperature power-law than possible using either data 
set alone. The radial temperature profiles of the disk models presented in Sections 3.2 and 
3.3 are shown in Figure [71 The temperature profile of the best fit power-law model is very 
similar to that of a classical, geometrically thin accretion/reprocessing disk (a = —0.75). 

As discussed earlier (Section 3.2), the accretion disk model requires an unrealistically 
large accretion rate, on the order of a few times 10~ 5 M yr -1 in order to fit the 
SED data. Simple estimates of the accretion rate for this star range from 5xl0~ 7 to 
0~ 6 yr" 1 based on mass loss data and upper limits to radio emission at 3.6 cm 



(IKurchakov et al. 20071 ; 



1 

Boehm fc Catala 19951 ; ISkinner et al. 19931 ). The necessity of an 



unrealistically large accretion rate in the accretion disk model is probably an artifact of the 
physically unrealistic assumption of a zero-thickness disk which causes the re-processing of 
radiation to be very inefficient. Hence, we do not favor an accretion disk model for the 
MWC 419 disk. 

A power law temperature profile with radius in which the radial exponent is a 
free parameter is a very good fit to the data, because such a simple model roughly 
accounts for the vertical scale height as a function of radius. The derived slope of -0.71 
i s consistent with a "flat" disk geom etry that one tends to find for Herbig Be stars 



flAcke et al. 2008 



Kraus et al. 2008h . 
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Radiol Temperature Distribution 



Radial distance (AU) 



Fig. 7. — The radial temperature profiles of various disk models. The dotted and dashed 
lines refer to the two cases of the classical-accretion disk model (Section 3.2; Table 1). The 
dotted-dashed line refers to the power-law temperature gradient disk model (Section 3.3; 
Table 1). 
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Earlier studies (IMonnier fc Millan-Gabet 20021 ; Eisner et al. 20041 ) using broadband 



usually single-wavelength, interferometric data, recognized a difference in the near-infrared 
size vs. luminosity behavior of high luminosity objects (pre-main sequence Be) compared 
to lower luminosity ones (pre-main sequence Ae), the former being more consistent with 
"classical disk" models. This has been revisited most recently by Vinkovic & Jurkic 
(2007), who use a model- independent comparison of visibility to scaled baseline and find 
a distinction between low-luminosity (< 10 3 L ) and high-luminosity (> 10 3 L ) YSO 
disks where the luminosity break point corresponds to approximate spectral type B3-B5. 
These authors then modelled this observable with a ring and halo model for low luminosity 
Herbigs, a halo alone for T Tauris and an accretion disk for high luminosity Herbigs. 

However, multi-wavelength interferometric studies have not always supported these 
conclusions when objects are modeled in detail. Kraus et al. (2008) observed the B6 star 
MWC 147 in H, K and N bands. They performed Monte Carlo modeling and found that 
the interferometric and SED date were not well fit with a standard, irradiated accretion 
disk alone, but were well fit with the standard disk plus emission from hot, optically thick 
gas within the innermost radius of the dust disk. In contrast, Acke et al. (2008) observed 
the B1.5 star MWC 297 in H, K and N, and found that they could not fit the data with 
a single accretion disk, even with the radial temperature exponent as a free parameter. 
Instead, they used a three component geometric model, with characteristic blackbody 
temperatures of 1700, 920 and 520 K. Our results show that MWC 419 (B8, 330 L ) has the 
disk characteristics of a high- luminosity object in the categories of Monnier & Millan-Gabet 
(2002) but a luminosity lower than the break point of 10 3 L Q identified by Vinkovic & Jurkic 
(2007). Examining Figure 2 of Vinkovic & Jurkic (2007) clearly shows MWC 419 fits well 
within the population of lower luminosity Herbigs in their model-independent comparison. 
Their physical interpretation of the low-luminosity Herbig group is an optically thick disk 
with an optically thin dust sublimation cavity and an optically thin dusty outflow. Our 
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multi-wavelength results do not support grouping MWC 419 in that physical model, but 
perhaps incorporating multi-wavelength data in the model-independent visibility groupings 
would have produced a different result. The new multi- wavelength data set presented here 
provides significant new information to aid in determining the physical conditions of these 
young stars, and our data show that the disk surrounding the B8 star MWC 419 is closer 
in physical characteristics to the more massive Be stars than to the Herbig Ae and T Tauri 
stars. 

Another result from the spectrally dispersed data is that the measured V 2 is marginally 
higher at 2.17/xm than at adjacent wavelengths. This V 2 feature is seen in all 4 independent 



data sets. As this is the wavelength of Br 7, and Brackett recombinati on emission 



mcs 



such as Br 7 and Br 10-20 are observed in the spectrum of MWC 419 (lHarvey 19841 ). we 
attribute the difference in line and continuum visibilities to a more compact Br 7 emitting 
region. Similar behavior in the apparent size in the Br 7 line compared to the neighboring 
con tinuum has been s e en before in spectr ally resolved interferometry from both VLTI and 



KI ( iKraus et al. 20081 : Eisner et al. 20091 ) . The derived uniform-disk diameter for the Br 7 



emitting region is 3.33±0.03 mas, which is ~ 4% less than the continuum size (3.46 mas) 
around Br 7 (Figure [2]). The significance level of this detection of the Br 7 emitting region 
with respect to the continuum region is ~ 2 a. The actual size of the Br 7 emission line 
region could be much smaller than the value reported here because of the coarse spectral 
resolution of our measurements. Detection of Br 7 emission inside the innermost dust radius 
suggests that the disk is optically thin in the inner region where atomic hydrogen gas exists. 



5. Summary 

This article reports the first milliarcsecond angular resolution observations of a 
HAeBe star (MWC 419) providing L-band, as well as simultaneously obtained K-band 
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data, both spectrally dispersed. This multi-wavelength observational capability is well 
suited to probing the temperature distribution in the inner regions of YSO disks, which 
is very important for distinguishing between models and gaining insight into the three 
dimensional geometry of the inner disk. Such measurements could distinguish discrete 
spatial distributions, such as dust-rims, from relatively-smooth spatial distributions, such 
as classical accretion disks, based on their distinct wavelength dependent disk sizes. In 
addition, interferometric measurements in the relatively unexplored L-band provide needed 
constraints to the disk/envelope geometry and temperature structure. 

Simple geometrical pole-on disk models are used to infer a linear relationship between 
the derived object size and wavelength in the 2-4 //m region, suggesting a simple physical 
model for the disk. The steep slopes of these linear relations imply that the disk is extended 
with a radial temperature gradient. We find that the accretion disk model of Hillenbrand 
(1992) derived from SED analysis does not fit our interferometric measurements. An 
updated accretion disk model with accretion rate 1.5 times larger and inner cavity 2.4 times 
larger fits the K band data well, but lies 3a above the L band data. However, both of these 
classical accretion disk models predict an unrealistically large accretion rate of ~ 3xlCT 5 
M yr _1 to fit the SED data. A power law temperature profile with a slightly shallower 
slope of -0.71 fits both the spectrally dispersed interferometric measurements and the SED 
satisfactorily, suggesting a relatively flat disk geometry for MWC 419. The measured disk 
size at Br 7 reveals the presence of compact emitting hydrogen gas in the inner regions 
of the disk. A more complete sample of YSO disk observations with adequate wavelength 
and (u,v) coverage, plus detailed radiative transfer modeling, are required to address the 
intriguing inner disk geometry in these sources 
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